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\ Abstract 

Mm 

CN \ The solving of the Schrodinger equation for a position-dependent mass quantum sys- 

tem is studied in two ways. First, it is found the interaction which must be apphed 
on a mass m{x) in order to supply it with a particular spectrum of energies. Second, 
given a specific potential V{x) acting on the mass m(x), the related spectrum is 
found. The method of solution is applied to a wide class of position-dependent mass 
Ch ' oscillators and the corresponding coherent states are constructed. The analytical 

^ . expressions of such position-dependent mass coherent states preserve the functional 

Q^l structure of the Glauber states. 

■ 1 Introduction 

■ The problem of calculating the energies of a quantum system endowed with position- 
^ . dependent mass m(x) and subjected to a given interaction represents an interface between 
O \ theoretical and applied physics. Its antecedent can be identified with the concept of 

effective mass, introduced in the forties to discuss the motion of electrons or holes in 
semiconductors [1]. Successfully applied in describing the formation of shallow energy 
levels due to impurities in crystals, the effective mass theory was strongly developed in 
the fifties [2]. Further insights were given in the calculation of superlattice band structures 
5^ \ for which the band edges and the masses are position dependent. In such context, it was 

stressed that the correct effective Hamiltonian consists of the kinetic term t-IP^, — ttI 

instead of the conventional expression |^ [3]. That is, the Hermiticity of the Hamiltonian 
is a part of the problem if the mass is not a constant. 

The subject has embraced potentials other than the periodic ones over the years. 
Indeed, the energy bands and periodic-like interactions appearing quite naturally in semi- 
conductor physics are substituted with point spectra and properly defined potentials in 
mathematical physics [4-13,21-28]. This new perspective has inspired intense activity in 
the looking for new exactly solvable potentials in Quantum Mechanics. Of particular in- 
terest, Susy-QM (factorization or intertwining method) [14-17] and group-theory [18,19] 
(see also [20]) represent the most fruitful approaches on the matter [21-28]. However, the 
literature focuses on mainly one of the two sides of the problem. Namely, in order to get a 
particular spectrum, the appropriate mass-function m{x) and potential V{x) are usually 
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looked for. A more deeper insight is necessary if one is interested on a particular mass 
m{x) subject to a previously defined interaction rather than looking for the recovering of 
a specific spectrum. 

In this work we analyze the two sides of the position-dependent mass problem by 
following the transformation scheme of the Schrodinger equation reported in [20] . In a first 
step the equation involving m{x) is mapped to the equation of a constant mass tjiq. After 
obtaining some general results, we study the eigenvalue problem connected with diverse 
position-dependent mass oscillators. In general, we distinguish between two fundamental 
kinds of oscillators. The first one is characterized by exhibiting the conventional set of 
equidistant energies hwo{n + 1/2), no matter the explicit form of m{x) or V{x). The 
oscillators of the second kind, on the other hand, are endowed with position-dependent 
mass m(x) and subjected to the conventional oscillator interaction Vosc{x) = rriQUjQX^ /2. 
The spectra of these last oscillators depend on the explicit form of the mass-function. In 
this way, we are able to compare the behavior of a quantum system of mass m{x) with 
that of a particle of mass mg when both of them are acted by the same oscillator-like 
potential. One of our motivations to analyze such oscillators is due to the fact that, as far 
as we know, there is a lack of results including the coherent states for position-dependent 
mass systems. 

Originally derived for electromagnetic fields [29], the features of the standard coher- 
ent states (Glauber states) are a consequence of the oscillator dynamical algebra [30]. 
They are usually constructed as eigenstates of the annihilation operator but are shown 
to minimize the uncertainty relation between position and momentum as well. A third 
property is that the Glauber states are displaced versions of the ground wavefunction. 
For other systems, generalized coherent states (CS) can be constructed through algebraic 
techniques (see e.g. [30,31]). In general, the CS show not all the three basic properties of 
the Glauber states. They have been recently studied in connection with non-linear Susy- 
algebras [32,33] (see also the reviews [34,35]), classical motion models for the Poschl- Teller 
potential [36], anharmonic vibrations in diatomic molecules [37], Landau levels [38] and 
the Penning trap [39]. With the present work we introduce some families of position- 
dependent mass oscillator coherent states. 

The paper is organized as follows. In Section [2] the Schrodinger equation of a position- 
dependent mass system is connected with the equation of a constant mass ttlq. The 
solutions are interrelated by a mapping for which the Hamiltonians are isospectral. The 
main difficulty is that the Hamiltonian of the mass mo includes an effective potential 
which, in general, makes the related equation as complicated to solve as the initial one. 
Here, the problem will be faced by either selecting the appropriate mass-function m{x) 
or by fixing the order in which m is entangled with P in the Hamiltonian. In Section [3] 
the previous general results are particularized to the harmonic oscillator potential. As 
a first result, it is shown that potentials behaving as a confining odd-root-law, In^, or 
sinh^ give rise to the quantum oscillator energies. The singular oscillator V = x'^ + otx~'^ 
is analyzed as a special case. On the other hand, it is also shown that the acting of an 
oscillator potential on m{x) involves the energies of a constant mass mo subject to either 
a confining even-power-law or the sinh^-like potentials. To deal with these last potentials 
it will be unavoidable the numerical approximation of the solutions. 
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In Section H] a position-dependent mass Hamiltonian is shown to be factorized by 
a couple of two mutually adjoint operators, the commutator of which depends on the 
explicit form of m{x). The initial Hamiltonian is then intertwined with a new one in 
such a way that they are isospectral. Factorization operators can be properly selected 
to work as ladders when acting on the eigenfunctions derived in Section [31 The related 
coherent states are constructed as eigenfunctions of the annihilation operator. These 
position-dependent mass CS are shown to have the same analytical form as the Glauber 
states. Moreover, they are also displaced versions of the ground state and minimize the 
uncertainty relation between P and X. Finally, in the very last section of the paper some 
concluding remarks are given. 



2 The eigenvalue equation 

Let us consider the one-dimensional Hamiltonian 

Ha = ^m^Pm^^Pm^ + V = Ka + V, 2a + 26 = -1 (1) 

where the mass m > and the potential V are functions of the position, Ka is the kinetic 
term of Ha and P fulfills [X, P] = ih, with X the position operator. We shall use to 
represent the domain of definition of Ha, i.e. V^""^ = Dom(Ha). 

In the position representation X = x and P = —ih-^, so we have: 

[fix),P]=^hnx), '^j-. (2) 
This last commutator allows us to express the Hamiltonian ([1]) as follows 

Ha =:^P' + ^(-)p-^ (4) - (2 + a){n^f] + V 

Zm 2m \ m J 2m Vm^/ -' ^^^^^ 



d? d 

= ai 

with 



ai-T^ + a2— + as. 



ai = -- — , a2 = -ai ( — ) = a' = ai ( ) \mm" - (2 + a)(m'f] + V. (4) 



( m' 

— , a2 = -ai — 
2m \ m 

In order to solve the eigenvalue equation 

Hai^ix) = Eij{x) (5) 
we first rewrite the functions ip to read 

^(x) = e3(^V(a;), (6) 
with g and ip two functions to be determined and such that 

\^lj{x)\^dx= [ \e'^''^ip{x)\^dx<+oo. (7) 
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Hence, from ([5]) one gets 



ai^" + {2aig' + a^W + {a^ [g" + [g'f] + a^g' + 03 - E}^ = 0. (8) 
Now, we introduce a change of the independent variable ruled by a bijection s as follows 

y = s{x), y X = s~'^{y). (9) 



The Jacobian of the transformation is given hj J = s'{x). If J 7^ at a point x, the 
inverse function theorem indicates that the map s is 1-1 and onto in some neighborhood 
of X (see e.g. [40], pp 91). In this way, to construct a well defined bijection s we first ask 
for the involved Jacobian s' to be free of zeros. On the other hand, let / be a function of 
X. Then we have: 

fix) = fis-\y)) = [f o s-'M ^ fM- (10) 

Thus, /* is the representation of the function / in the y-space. In a similar manner 
we find that / is the representation of in the a;-space: / = [/^, o s]. Hereafter, and 
whenever there be no confusion, we drop the subindex "*" from the functions in the 
^-representation. The straightforward calculation departing from Eq. ([8]) leads to 

ai(i/')V + M' + (2«i^7' + «2)2/']0 + {aiW + [gT] + ^29' + as - Ejcp = 0, (11) 

where / = df/dy. This last equation acquires a simple form if the coefficients of if and 
if are respectively a constant Cg (expressed in appropriate units) and zero. Thereby one 
has a system of equations 

a,{yT = cl, a,y" + i2a^g' + a2)y' = 0, (12) 

the solution of which defines the form of g and y in terms of the mass position-dependence: 

1/4 



g{x) = In 



mix) 



mo 



y 



(13) 



Here mg and yo are integration constants (we shall take, for simplicity, yo = 0) and the 
constant cq, introduced in f|T2l) . has been written as Cq = ih/ ^2mo- These last results in 
f fTT]) reduce the initial eigenvalue equation ([5]) to the following one: 



where the function 



2m(] 



^{y) = E^{y) 



(a) 



eff 



V 



2m^ 



— ha mm 
4 



/\2 



{m) 



(14) 



(15) 



plays the role of an "effective potential" and depends on the explicit expressions for the 
mass m and the initial potential V , both of them in the ^-representation. In particular, 
if the mass function m is a constant then we have V^^'' = V . In general, m{x) could lead 
to a very complex function V^^ {y) for which the new equation is as complicated to 
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solve as the initial one (Q. Hence, at this stage, the main simplification is the avoiding of 
undesirable mass factors in the derivative term so that techniques to solve the conventional 
Schrodinger equation can be applied. 

Given a solution ip^ of (fT4l) . according with (j6]), and (fT3|) . the function ip is 



mo 



1/2 



(16) 



Therefore we have 



x)\ dx = / J\(p^{s{x))\ dx < — > / \ip^{y)\ dy, (17) 



(a) 



with = Dom{H^^). That is, by getting the square- integrable eigenfunctions of 
one is able to obtain the square-integrable eigenfunctions of Ha and viceversa. Moreover, 
from dS]) and f[T^ we notice that (p^, and if) share the same eigenvalue E. Thus, Ha and 
iJgg are isospectral operators; we write Sp(//a) = Sp(i^gg ). 

Notice that equations (ITS!) and (fT6|) are consistent with the results reported in [5,23,25]. 
With regard to our approach, there is yet a couple of special cases leading to further 
simplifications. Namely, one can get V^^{y) = V{y) by selecting the appropriate function 
m{x) or by properly fixing the value of a, as we are going to show. 

2.1 Mass-dependent null terms (MDNT) 

Let us look for a mass function m such that V^^^'' — \^ = in Eq. f[T^ . Thus, we should 
solve the non-linear, second order differential equation: 

cimm" + C2{m!Y = 0, 2ci = — C2 = — h 2a. (18) 

16 2 

A brief examination yields 

m{x; a) = mo{xo + Ax)-^/(^+^'^), a ^ -3/4 (19) 

as the simplest solution with xq and A constants to be fixed. We have to distinguish 
between two general cases: 

I) If a < —3/4 then m{x; a) has a zero at x = to = —Xq/X 
II) If a > —3/4 then m(a;; a) is singular at x = to 



The first case will be omitted to avoid ill defined operators Ha and unappropriate map- 
pings s as well. Indeed, if a < —3/4, the kinetic term Ka in Eq. ([T]) diverges and the 
Jacobian J in (|T6l) is zero at x = to- On the other hand, for a > —3/4 the integrability 
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of ^m{x] a) in Eq. (fT3|) depends on the value of a. In particular, ifa = ao = — 1/4 then 
the mapping x ^ y^o) is ruled by the function 

ln(xo + Ax) 

y{0) = S(o){x) = , x>to, (20) 

with 

m,o,{x)^mix;ao) = j^-^. (21) 

This last expression of m{x) is connected with the revival wave-packets in a position- 
dependent mass infinite well [10]. Here, the Jacobian reads J(o) = l/(xo + Ax), so that 
the bijection S(o) is well defined for all x > to and arbitrary real values of xq and A. Then, 
in general p("o) C [to, +oo) and ^^g^"'' C M. The explicit form of the domains of definition 
and "D^g ^ depend on V while the inverse function reads 

x = s^^,{y) = '-^^^. (22) 

On the other hand, if a 7^ Oq the new variable is given by 

(xo + Ax)^ f l + Aa \ 

y = s{x;a) = ^-^, V=[j^^). (23) 

The appropriate mapping s(x; a) is fixed by looking for the values of a such that either 
1] = 2n + 1 or 1]'^ = 2n + 1, n = 1,2,... In the former case we arrive at the discrete 
set integrated by the points — = "~l)~8'~i'''' ^ ~4' However, each one of these 
possible values of a belongs to the case (I) discussed above and must be omitted. Now, 
if ri~^ = 2?T, -|- 1 we obtain the points = 0, — |, — |, • ■ ■ > —j, which belong to the case 
(II) we are interested on. Henceforth, the mapping x is established from Eq. fl23l) 

as follows 

?/(„) = s(x; a„) = S(„)(x) = (^—y—^ (^0 + Ax)^/^^"+^\ 7^ ao, (24) 

while the corresponding inverse transformation is ruled by 

1 



A 



\ \ 2n+l 

AZ/(n) 



2n + 1 



Xo 



an 7^ ao- (25) 



The expression for the mass-function (fT9l) in terms of a„ 7^ ao reduces to 



(xq + Ax)4"/(2-+l) ' 



m(„)(x) = m(x; a„) = — — ^_\!^r,/('?.r,-l-^^ ^ n G N. (26) 



Remark that J(„) = (xo + Ax) 2n/(2n+i)_ Hence J(n) 7^ for all x G M and arbitrary real 

i7 



values of xq and A. As a consequence C M and P^g"^ C M. To embrace a„>o and ao 



into the same notation let us introduce the set 

1 

4n 



1 — 71 

A= {ao = -1/A,an = — — \, nGN. (27) 



Then, if a & A the position dependent mass operator Ha„ is mapped to a conventional 
Hamiltonian in the ?/(ri)-representation and viceversa (see Tabled]): 

/f->^-(^)^ + n«„,). (28) 

with Sp(ifa„) = Sp (if We shall take full advantage of this last property in the next 
sections. 



2.2 Mass-independent null terms (MINT) 

A simple inspection to equation f[T^ shows that v}f^^^^\y) = y{y), no matter the explicit 
form of the mass function m{x) -assuming this last is well defined-. That is, by fixing 
a = —1/4 we get: 

In particular, if m = m(o) then H in ( l29i) is the same as H'^°'°^ with = Dom(if_i/4) = 

p(ao) g^j^^ Dom(iy) = T>'^^\ A similar situation occurs if m = m(„) (see Tabled]). 

Besides the mass-functions derived in the previous section, a regular expression for m 
has been recently introduced in [27] (see also [7]). This is given by the function rriR. 

1 + [Xxy 

with 

, , arcsinh(Aa;) i, , sinh(Av) 
sr{x) = s~^\y) = j-^. (31) 



The corresponding Jacobian J = l/yl + (Ax)^ is nonzero for all x G M and arbitrary 
values of A. Hence P^"^/^) C MandP^g^^^^ C M. The main aspects of these last results are 
summarized in Table dl As a final remark, the mass (130]) appeared in the construction of 
the Wigner functions connected with a class of position-dependent oscillators [13]. Other 
interesting mass-functions are 

W + X 



m^ix) = y ^j^^2 j ' ^cix) = cx . (32) 

They have been already studied in [22] and recently in e.g. [6]. Notice that Jc{x = 0) = 
and Ju, 7^ Vx G M and w > 0. In the next sections we shall study specific forms of 
the potentials V{x) and V^..{y) which represent oscillator-like interactions for a position- 
dependent mass quantum system. 
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\2m.oJ rf^ 



ao 



1 -l/4p -l/2p -1/4 
2*^(0) -^"^(0) -^"^(0) 



Pf-^o) C [to,+oo) 



ai = 

m(i)(x) = ^ 



Xy(i) ^ HxQ + Xx^/^ 



7-)(ai) ^ 



m(„)(x) -see Eq. (^5]) - 
j/(„) -see Eg. ([HD- 



1 On 



p m^^-r^"" p mf::, + v{x) 



2mo y 



^(y(n)) 



m{x) > 

2/ = J(m/moy^^dx 



H-i/4 = \ m~^/^Pm-^/^Pm-^/^ + V{3 
p(-i/4) c R 



H 



2mo y dy 



en — 



V{y) 



Table 1: Special mass functions m{x) and orderings of the kinetic term Ka{x) leading to the Hamil- 
tonians Ha 
fixed by V. 



tonians Ha and h'^^ , with V^s\y) = ^iv) and Sp(iJa) = ^H^^). In aU cases the definite domain is 



3 Two kinds of position-dependent mass oscillators 

We are going to work with the eigenvalue equation (fT^ such that V^^ =V hy either the 
MDNT or the MINT cases described in the previous sections. Although our approach 
holds for any well defined potential V , we shall focus on the linear harmonic oscillator in 
two general situations: 



i) Departing from a given interaction V{x) and a mass function m{x) we arrive at the 
conventional linear harmonic oscillator problem in the ^/-representation. That is, the 

2 

new potential reads = — y-^?/^, with ujq the natural frequency of oscillation. 

Since V{x) and are isospectral they share the eigenvalues defined by En = 

huJo{n + 1/2), n = 0, 1, 2, . . . We shall refer to these potentials as oscillators of the 
first kind. 

2 

ii) Departing from the linear harmonic oscillator interaction V{x) = "^"^^ and a 
mass function m{x) we arrive at the eigenvalue equation connected with the new 
potential V^{y). Since V{x) and V*{y) are isospectral we solve the (conventional) 
Schrodinger equation in the ^-representation to construct the solutions of the initial 
oscillator-like, position-dependent mass problem. We shall refer to these potentials 
as oscillators of the second kind. 



3.1 Oscillators of the first kind 

2 

Let us take K<(y) = ^^^^y^ as the ^-representation of the initial potential V{x). Then 
= M and all the mappings MDNT and MINT can be applied (see Table [I]). It is 
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convenient to introduce a dimensionless notation as follows: 



2 rfy2 2 



Then, the solutions read 



<^(y) = o, J = y[ — ] =y(^^ e 



E 

hjJn 



^nij) = — -, -n.„ y = — 1 e-'^ ' - — e ' En = n + -. 



(33) 



(34) 



Next, we are going to solve the initial position-dependent mass problem in terms of these 
results. 



3.1.1 MDNT case 



Let m{x;a) be the mass function with a & A, that is m = m(„), n 
equations fl20l) and flMl) we know that the initial potential reads 

^^^2 [ (2n + l)2(xo + Ax)(5^, nGN, XG 



0,1,2,... 



ln^(xo + Ax), 



n 



0, X G [to, +00) 



From 



(35) 



A dimensional analysis shows that A = Aqo;, with Aq a constant (Aq = 1 for simplicity). 
The behavior of potential ( |35i) for = 1 and n = 2 is contrasted with the well known curve 
of the harmonic oscillator potential in Fig.[T^; the case n = is depicted in Fig. [lb. In both 
cases, as we have previously noted, the involved spectrum is given by E^ = hwo^k +1/2), 
= 0, 1, 2, . . ., while their eigenf unctions respectively read 



Hk[{2n + l){xo + ax] 



(2n + l) 



(xo + ax) a/2%V2^ 



exp 



(2n + l) 



-(xo + ax)(2"+i) 



.n eN, 



and 



Hk[\n{xo + ax)] 



n = 0, 



(36) 



(37) 



with X running in the domains indicated in fl35l) . Thus, the energy spectrum of a position- 
dependent mass quantum system which is subject to either the action of a confining odd- 
root-law potential V(„)(x) oc (ax)'^^^'^"'^^\ n G N, or to a square- logarithmic interaction 
V(o)(x) oc In^(ax), is ruled by the quantization of the conventional harmonic oscillator 
energy if the mass- function is respectively taken as m(„) or m(o). 



3.1.2 MINT case 

Let mFi{x) be the mass function with A = a. Then, the spectrum of the potential 



v{x) = [v;os](x) 



fkJn 



arcsinh (ax) 



(38) 



is given by Ek = huJo{k + 1/2), k = 0,1,2, .. ., and the involved eigenfunctions read 

1/4 



mix] 



_22%(A;!)2mo 



Hk[aTcsmh{ax)] e'^ arcsinh {ax)^ k = 0,l,2,... (39) 



In Fig. [T] the global behavior of potential (l38l) is shown in contrast with the curve of the 
harmonic oscillator one. 




(a) 



Osc 
ASh 
n=l 
n=2 



2 . 5 



1.5- 



. 5 




-10 1 



(b) 



Figure 1: (a) The odd root-law potential ((35)) with n = 1, n = 2 and xq = besides the regular one 
(ASh) defined in Eq. ([55)) . The harmonic oscillator potential (Osc) is depicted as a reference, (b) The 
square-logarithmic potential (Log) defined in Eq. ()35p with xq — I and — [— l,-|-oo). Potentials 

(Ash) and (n=l) as well as the first three energy levels (E) are also depicted. Vertical and horizontal 
axis are respectively in Hluq and dimensionless units. 



As we can see, one is able to identify the kind of interaction V{x) which has to be 
applied to a quantum system of position-dependent mass m{x) in order supply it with a 
specific, well known, spectrum Sp(//*^"^). For instance, if Sp(if*^"^) = {huQ^n + 1/2)}^^, 
we have shown that the system has to be subject to potentials behaving as a confining 
odd-root-law. In or sinh , whenever the mass-function is respectively defined by ( I2T1) . ( !26l) 
or ( l30l) . A more deeper insight is necessary if one is interested on a position-dependent 
mass m{x), subject to a particular interaction V{x), rather than in the recovering of 
a given spectrum. That is, what the sort of the spectrum is expected by applying an 
oscillator-like interaction to a quantum system of mass- function m{x)7 The problem is 
going to be faced in the next section. 



3.2 Oscillators of the second kind 



In this section we analyze the effects on the energy spectrum produced by a position 
dependence of the mass. In other words, how different is the spectrum of a system of 
mass m{x) from that of a particle of mass mo when both of them are under the action 
of the same potential K(x)? As before, we shall focus on the simplest case of the linear 
harmonic oscillator interaction. 



Let V;< 



[X 



2 



be the initial potential. Notice that i5om(V^f 



requires 



= M. However C [tQ,+oo), so that a must be different from qq (see Tabled]). 

The case a = will be analyzed in Section 13.2.31 
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3.2.1 MDNT case 



Let m 



in) 



be the mass function with n G N and X = a. The initial potential Vq, 



behaves as an even-power-law function in the ?/(„)-space (see Fig. [2]): 



V;(i/;n) = [K. 



o s 



in) 



ay 

2n + 1 



2n+l 



- Xo 



n eN 



(40) 



where the label "(n)" has been dropped from the y-coordinate for simplicity. Hereafter 
we shall take xq = 0. Notice that K(l/; n) — as n — +oo and K(l/; n) — hw^^ayY /2 
as n — > 0. Thus, the family of potentials (HOl) is delimited by the free particle and the 
harmonic oscillator potentials (remember that n = and n — >■ +oo are forbidden in 
Eq. [21]). Such a behavior is shown in Figure [2J 
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Figure 2: Three members of the family of power-law potentials (l40l) . The conventional harmonic 
oscillator potential (Osc.) is recovered for the forbidden value n = and the family goes to the free 
particle case for n +00. 

Let us emphasize that, although this kind of potentials is not analytically solvable, 
they have deserved special attention in pedagogical as well as in research papers over the 
years. For instance, their WKB energy levels have been shown to depend on the power of 
the potential [41] and the involved scale invariance has been studied in terms of the Lie 
method [42] . The polarizability of a particle in a power-law potential due to the presence 
of a constant force and the wave packet revivals in such potentials, on the other hand, 
have been exhaustively studied in [43] and [44] respectively. 

In general, the roots of V^{y]n) = Ksc(2/) define a region {—yc,yc) C ^^es"^ which 
the potential K grows up slower than Vosc- The geometry of these last curves in {—yc, yc) 
suggest the spectrum of will be integrated by energy levels which are below the cor- 
responding oscillator energies. The behavior of the curves in the complementary region 
\ (~2/c, Vc) is such that the energy levels are expected to be above the oscillator ones. 
To verify our statement let us calculate the eigenvalues of K(2/; n) by means of the energy 
quantization condition of the WKB method: 

/+yo 
v/2mo[^(n) - V{y- n)] dy = 7ih{k + 1/2), A; = 0,1, 2,... (41) 
-yo 

with ±?/o = ±(^^^)(2E(n))^/*-^"'"'"^^ the classical (symmetric) turning points and E{n) the 
energy connected with the potential V^{y]n) for a given n G N. The change of variable 
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y = hqz reduces the integral equation (HTl) to (compare with [41] and [42]): 



Ek{n) 



hujn 



71 (fc+ 1/2) 

,Jn (2n + 1) 



2n + l 
n + 1 



where the constant 



Jr. 



2(^ + i)r(fS) 



(42) 



(43) 



is such that j„ — > 2 as n +oo and jo = 7r/2. Figure [3] shows the spectrum curves 
of three members of the family fHOj) compared with the energy spectrum curve of the 
harmonic oscillator. Notice that the energy levels become closer to each other as the label 
n increases (free particle case). That is, if n » 1 then Ek{n) oc [{k + l/2)/n]^. On 
the other hand, for the forbidden value n = we have the oscillator spectrum Ek{0) = 
huo{k + 1/2), as it was expected. The corresponding set of eigenf unctions, in turn, can 
be numerically constructed or analyzed by using improved versions of the WKB method 
hke that discussed in [45]. 
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Figure 3: (Left) The spectrum curves Ek{n) of the power-law potential (|40|) for n = 1,2, 3, besides the 
spectrum curve of the harmonic oscillator (Osc); all of them are depicted in huo units as a function of 
k. Notice the points in which Ek{n) = Ek{Q)- (Right) Details of the first four energy levels. 

Now, let us take one of the allowed values of n. The root of equation Ek{n) = Ek{0) 
is given by 

kc{n) = i 

The ceiling function \kc] of kc{n) defines a subset of Sp(K(|/; n)) = {Ek{n)}^^ for which 
Ek{n) < -Efc(O) Vfc < \kc~] ■ The larger the value of n the bigger the set of eigenvalues E^iji) 
dehmited by -Efc(O). The complementary set is then such that Ek{n) > Ek{0) Vfc > \kc]; 
details are shown in Figure El 

In conclusion, the oscillator of the second kind defined by the pair (Ksc, "^(n)), n eN, 
shares its spectrum with a particle of mass mo subject to an even-power-law potential of 
the form K(l/(n); n) oc [?/(„)/(2n + 1)]*""*"^. When contrasted with a conventional oscillator 
of mass mo, the energy spectrum of the pair (Ksc,'"^(n)) is a distorted version of the 
oscillator one. The shape and amount of the distortion are respectively dictated by Eq. 
(Il2|l and kdn), as this last was defined in (14^ . That is, the distortion is stronger for 
larger values of | [fed — k\ 'va. equation 



(2n+ l)r( 



1 



4n.+2 ' 



n+l - 
2n+l' 



Zn + l 



(44) 
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4 
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5 


10.08292 


10.11165 


6 


12.63890 


12.66657 


7 


15.37683 


15.40365 


8 


18.28821 


18.31431 


9 


21.36592 


21.39141 



Table 2: The first ten energy levels of the potential sinh^(y) calculated numerically from the WKB 
transcendental equation (|47l) and directly from the Schrodinger equation (|46p . 



3.2.2 MINT case 

Let m^{x) be the mass function with A = a. The potential in the ^/-representation reads 



K(y) = [Ksc0.^1](y) 



■ sinh [ay). 



Here, the (dimensionless) Schrodinger equation to solve is 



df 



+ sinh y 



(45) 



(46) 



As in the previous case, the energy quantization condition (HT!) gives an accurate ap- 
proximation to the eigenvalues E of the energy. With the classical turning points ±?/o = 
± arcsinh (\/2E) /a, one arrives at the following transcendental equation: 



2E Fe z arcsinh V 2E 



1 

2E 



vr 



^-(A; + l/2) 



(47) 



where 



m) 



{l-m^An^ef^de 



is the Elliptic Integral of the Second Kind (see e.g. [46]). The roots E^ of fHTI) can be 
evaluated numerically by using conventional algorithms. In Table [2] we show some of the 
first values of E^ compared with those obtained from a direct, numerical integration of the 
Schrodinger equation (H6!) . The corresponding probability densities |(/9^,(y)p are plotted 
in Figure m contrasted with their partners |'?/'(x)p in the ^-representation. 

In this case the geometry of the curves K(|/) and Vosdy) is in correspondence with 
the fact that all the energy eigenvalues of K are above the related energy levels of Vosc- 
Indeed, around the origin one has K > Ksc, so that Eq > 0.5/ia;o, as expected (see 
Table [1]). For an arbitrary excited level Ek, the distortion is as strong as fast is the 



13 



Figure 4: (Left) The potential sinh^(y) and the probabihty densities of its three first wavefunctions 
together with the corresponding energy levels. The oscillator potential curve (dashed) is included as 
a reference. (Right) Probability densities of the first three wavefunctions of a second kind, position- 
dependent mass oscillator which shares its spectrum with the sinh^(y) potential. In both cases the 
vertical and horizontal axis are respectively in Hujq and dimensionless units. 

growing up of K — Vosc- In conclusion, the oscillator of the second kind (Ksc,''^_r) shares 
its spectrum with a particle of mass mo subject to the sinh^ potential (see Fig. Hj). The 
spectrum, in turn, is a strong distorted version of the conventional oscillator's one. 



3.2.3 The squeezed oscillator 

Let us consider the potential 



Xf) + ax 



— {xq + ax) 



+ 2(1 



X > t 







(4J 



with a defined in ( l33l) . Xq a dimensionless constant and Dom(V^q) = [to,+C)o). This 
potential is often refered as the "singular oscillator" because its singularity at x = to- The 
conventional expression V^^^x) = mo^o {y + j , with g in units of the square of distance 

and shifted by — is recovered from ( l48l) with xo + ax = z and a^g = \/2/A. Here, 

we prefer to call it sqeezed oscillator because its domain of definition is the result of a 
L.H.S. 'squeezing' of R in terms of s~^, as it was established in the previous sections. If the 
mass function m{x) is a constant mo the involved (dimensionless) Schrodinger equation 



(l-.)^2(v/2-l) 



V9 = 2Ev9, 



Xq + ax 



(49) 



can be solved in terms of confluent hypergeometric functions by means of the appropriate 
transformation (see e.g. [20]). Indeed, the mapping 
the following Kummer equation 



z e 



-z2/4 



uiz), z 



2x, leads to 



V2 



X 



dx 



-u 



E U = 0, 



(50) 
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Figure 5: The squeezed oscillator ([^S]) and its first three probability densities. Vertical and horizontal 
axis are in hujQ and dimensionless units respectively. 



with . 
read 



'l+\^)/2. Thereby, the physical solutions for En = }Tioo{n+l/2), = 0, 1, 2, . . ., 



Cn 2; 2 e 



1/2 

z 2 e 



(51) 



2i/v^r(n+l+l/V2)^ 

with L'^n\x) the Generalized Laguerre Polynomials [46]. So that, for a const ant- mass 
quantum system, the one- dimensional potential ( HHi) shares its spectrum with the con- 
ventional linear harmonic oscillator. Moreover, it is well known that formulae (15T|) can 
be also algebraically obtained in terms of the sm(1, 1) Lie algebra (see e.g. [30], pp 217). 
The squeezed oscillator is shown in Figure [5l together with some of the corresponding 
probability densities. 

If m{x) is not a constant then the energy spectrum of the quantum system is modified, 
as we have previously verified. According to Table [H any of the masses f l2T|) . (l26l) or (!30!) 
allows the mapping to the y-space. First let us consider the case m = m(o). The potential 
IHll is mapped to the following one 



K(l/(0)) = [V;qOS(o5](Z/(0)) 



sinh^(a?/(o)) + 



1 - V2 



(52) 



which, up to an additive constant, is the same as the potential reported in Eq. (H5l) . 
Thereby, we have shown that a quantum system endowed with mass m(o)(a;) and acted 
by the oscillator-like potential (HHI) shares its spectrum with a particle of mass mo which 
is under the action of the potential sinh (y(o))- In comparison with a constant mass 
quantum oscillator, we realize that the presence of m(o)(x) distorts the ground energy 
level of the oscillator-like system from 0.5 to ~ 0.6 energy units {hwo), the second one 
from 1.5 to ~ 1.9 and so on. The higher the level of excitation of the system the stronger 
the distortion of the spectrum. 

Notice that the system studied in Section 13.2.21 behaves in a similar manner, so there 
exists a clear relationship between position-dependent mass systems: different masses 
combined with appropriate interactions give rise to the same spectrum. In this case, the 



oscillators of the second kind defined by the pairs (Ksc, f^R) and (Kq, 'm(o)) are isospectral 
(see Table |2]). On the other hand, we have another pair of oscillators of the second 
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kind (V^q,m(„)) and (Ksq, ""^ij), which are respectively isospectral with the const ant- mass 
potentials 



tujJn 



n)i 



2n + l 

ay{n) 



2n+l 



(^y[n) 

2n + l 



2n+l 



+ 2(1 - V2) > , n e N (53) 



and 



huJn 



1 



1 2 



xq + sinh ay 



— (xo + sinh ay) 



2(1 - V2) 



(54) 



Each one of these last potentials shows a spectrum which is a distorted version of Sp(V^q) = 
Sp(Vosc)- In summary, given an interaction represented by V{x), the spectrum of a 
position-dependent mass quantum system is a distorted version of the spectrum of a 
particle of mass ttlq subject to the same interaction. As we have realized, the degree of 
distortion depends directly on the explicit position dependence of the involved mass. 



4 Factorization and Coherent States 

Once we have constructed the solvable position-dependent mass Hamiltonians Ha, one can 
look for the appropriate factorization operators. The presence of m{x) in Ha makes nec- 
essary a refinement of the factorization (see e.g. [20] and [47]). As usual, the factorization 
operators intertwine the initial Hamiltonian with a set of new exactly solvable energy-like 
operators Ha [15, 16]. However, in general they do not act as ladder operators on the 
eigenfunctions of neither Ha nor Ha- In the case of position-dependent mass oscillators 
of the first kind, the factorization operators act in a ladder form if their commutator is 
the appropriate constant. Then, as we are going to show, it can be constructed a set of 
position-dependent mass coherent states. 



4.1 The position-dependent mass factorization 

Let A and B be the following operators 

A = — ^-^m'^Pm^ + (3, B = -^m'Tm" + f3, A^ = B, (55) 
\/2 V2 

with (3 a function of the position operator X. We want to work with A and B as the 
factorization operators of Ha. In this regard, it is important to stress that most of the 
literature pay attention on a specific ordering of m and P. Namely, it is usual to take 
a = and 6 = —1/2 so that the kinetic part of Ha reads |P^P, with the corresponding 
simplification of A and B (see e.g. [21,22,28]). Here, we shall use the operators (155!) 
with no a priori assumption on the ordering of m and P. In this way, the results already 
reported will be included as particular cases. 

If A and B factorize the Hamiltonian ([T]) in a refined way [47] then one has 

Ha = AB + e, (56) 
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and (3 fulfills a Riccati equation in the position representation: 



V -e 



h 



'2m 



2 a + 



4 



m 



m 



(57) 



where e is a constant (in energy units) to be fixed. For arbitrary m and (3 the product 
between the factorization operators obeys the commutation rule: 



[A,B] 



m 



3 V 4 



3(m 



/\2 



mm 



2h 
/2m 



Therefore we have a new operator H^, defined as follows 

Ha = BA + e = Ka + V, V:=V-[A,B], 
which is intertwined with Ha by means of the factorization operators: 

H„B = BH„. H„A = AH„. 



(5J 



(59) 



(60) 



The relevance of these last relationships is clear by noticing that, if is an eigenfunction 
of Ha with eigenvalue E (see Eq. [5]), then ifj oc Bip 7^ solves the new eigenvalue equation 



Ha^l){x) = Eij{x) 



(61) 



Moreover, it is easy to verify that a normalized wavefunction leads to IV'P oc E — e. 
Then, the new set {ip = Bip/(E — e)^/^ \ E ^ e} consists of normalized eigenfunctions 
of Ha belonging to the eigenvalues {E} = Sp{Ha). Now, let ipt be a function which is 
orthogonal to the set {ip}, i.e., {4',i'e) oc (V^, Aip^) = 0. Since 7^ we have Aip^ = and 
necessarily Haipe = e^e- The involved solution reads 



(62) 



with Ce a constant of integration. If {ip^,ip^) < 00 then Sp{Ha) = Sp{Ha) U {e}. 

The previous derivations considered e ^ Sp{Ha). To include the case e G Sp(ifa) let 
us assume that the solution of BipM = 0, given by 



(63) 



is a square-integrable function. In this way ^/'m is the wavefunction of Ha belonging to 
the eigenvalue E = e. As a consequence, there is no element in {^p} constructed from ipM 
via the relationships ( !60l) . The corresponding function ipM must be obtained as a solution 
of BAipM = (see Eqs. [591 and [6T]) . There are two possible cases: 

1) If AipM = 0, then ipM has the same form as the function defined in (l62l) . However, 
ii ipM € L'^iV^''^), from Eq. fl63l) one notices that ipM oc m^^'^/ipM is not square-integrable. 
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2) If AipM 7^ and B{Ai/jm) = 0, one can take AipM = i^M such that B{iIjm) = 0. 
Then, because ip and ipM are orthogonal, we have {%Ij,iPm) = (V')^V'm) = ii^^i^Ai) = 0. 
Thus, ip and ^/'Af are orthogonal and Sp(ifa) = Sp(ifa), with e G Sp(ifa). 

At this stage it is important to stress that Sp{Ha) and {ip} can be obtained by means of 
the transformations introduced in Section [2l Thereby, we get a wide family of isospectral 
operators if, for instance, H^"""^ is the MDNT-Hamiltonian defined in ( !28l) . That is, we 
have : 

SpiHaJ = Sp(if.J ^ Sp(/f("")) = Sp(/J("")) (64) 

where H^"""^ is the Hamiltonian intertwined with in the ?/(„)-representation. The 

same can be said about the MINT-Hamiltonian H defined in ( l29l) . In this context, it will 
be profitable to decompose the commutator ( l58l) in the MDNT and MINT cases: 



h / 1 

^3/2 I 4 



1 2 



m 



— 

m 



21)2 

/?' (MDNT) 

m 



(MINT) 



(65) 



4.2 Position-dependent mass ladder operators 



Let us consider the Hamiltonian of an oscillator of the first kind Ha- In advance we know 
that Sp(ifa) = {E„ = hujQ^n + l/2)}+^o, whichever the MDNT or the MINT case we are 
dealing with (see Section [3]). To get the simplest form for the corresponding annihilation 
and creation operators let us take [A, B] = —hwQ. Then we have Ha = Ha + hujQ. That 
is. Ha differs from Ha only in the zero of the potential. This physical equivalence and the 
intertwining relationships (160|) make clear the roles played by the factorization operators: 



A{Ha + hwo) = HaA, B{Ha - fuJo) ■■ 

Therefore Aip^ oc ipn+i and Bip^ oc ipn-i if Haipn = E^ip^ 
[A, B] = —hujQ in fl58|) leads to the following /3-function 



HaB. (66) 
Now, the substitution of 



72 



^ h 



V2 



m 



m 



3/2 



(67) 



Here /5o is an integration constant which will be omitted in the sequel. The identification 
e = Itwo/2, after introducing fl67|) in the Riccati equation fl571) . allows to write the potential 
V in terms of the /9-function: 



2h 



V2 



a + 



h 



2V2 



// 



mm 



im 



i\2 



(6J 



The straightforward calculation shows that this last expression is reduced in both the 
MDNT and the MINT cases to the same simple form: 



2 



m 



(69) 
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Remark that this expression for the potential is consistent with our transformations in 

2 

Sections [2] and [31 Indeed, since K {y) = y"^ has been given as the initial potential, its 
^-representation reads 

2 

V{x) = [K o s]{x) = V;(s(x)) = ^(s(a;))2 (70) 

with y = s{x) given in (fT3l) . In summary, we have shown that A and B, as they are defined 
in fl55ti57p . are nothing but creation and annihilation operators if their commutator (l58l) 
is constrained to be a constant equal to the separation between the energy levels of Ha. 
The same condition allows to identify a quadratic potential V, expressed in terms of 
the mass-function, which is consistent with the transformations defined in the previous 
sections. 



4.3 Position-dependent mass coherent states 



To take full advantage of the results derived in the previous sections let us rewrite the 
factorization operators as follows (compare with [21,22,28]): 



A 



B 



'2m 

h 

Jim 



d 

dx 
d 

dx 



(Inm)' 
4 

(Inm)' 



+ 



^0 

V2 



m 



+ / m^^'^dr 
V2. 



(71) 
(72) 



where we have used ( l55l) and ( 1671) . The operator B in the y-space is then given by 



h d 



+ 



1/2 



h 



y/2m^dy \ 2 J 



V32mo V dy 



d In m* 



(73) 



and a similar expression for A^, obtained from (173!) by changing the sign of the first and 
third terms. Finally, in the dimensionless notation of Eq. (!33|) we get 



d 

dy 



+ y 



d 

dy 



In 



1/2 



(74) 



and A^, = A*a/^o/2. Then, from fl55|) we have 

H,, = {2/huJo)Ha, = A,B* + 1, [A„ B,] = -2 (75) 

The action of 11^^, , A^, and B^, on ip in the dimensionless ^-representation is as follows: 

Ha*^/'* = Jy^(-V3 + yV), A*?/;^ = jy^a+v9^, B^V^^ = J^^a^v^^. (76) 

The Jacobian J is defined in Eq. (fT6l) and a_ (a+) is the conventional annihilation 
(creation) operator of the linear oscillator in the y-representation 



A. 

dy 



+ y, 



2N 



(77) 
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with the Fock's number operator. Hereafter we shall omit the "*-notation" . 

In order to construct a set of coherent states as eigenfunctions of B we first take an 
arbitrary linear combination B of the wavefunctions ipn associated with Hq: 

oo 

e = J2ckiJk. (78) 
The action of B on this last function reads 



B& = J^/^J2ckV2kipk-i. (79) 

k=0 

We look for the functions B fulfilling BB = zQ, z G C. The straightforward calculation 
leads to a recurrence relation which is satisfied by the coefficients c^. The root is found 
to be Cfc = z''cQ/Vk&. As usual, the coefficient Cq is fixed by the normalization of B and 
we finally arrive at the familiar expression: 

1 2 1 ^ ^ 2^ 1 \z\ ^'^-\- 1 \^\ — 1 

where z stands for the complex conjugation of z and we have used the Baker-Campbell- 
Hausdorff formula e^e^ = exp(A + B + ^[A, B]) to recover the displacement operator 



D{z) := = e^"+/2e-^'^-/2e-l^l (81) 

Thereby, since 6z{y) '■= D{z)ipQ{y) is a conventional constant-mass coherent state in 
the y-space we conclude that its partner in the x-coordinates B2 = J^^'^Oz is a position- 
dependent mass coherent state, defined in terms of the annihilation operator B. Explicitly, 
Qz{,x) is given by 

B.(x) = ("lI^M^^ ' f2 Ms{x)). (82) 

This last result is the general expression for the coherent states of any of the oscillators 
of the first kind introduced in Section [3?T1 In this context, let one of these oscillators be 
in the state B^. The probability of getting E„ = 2n + 1 as the result of a measurement of 
the energy is ruled by the Poisson distribution: 

Pn(e2) = |(^„,B2)p = LL.e-NV2. (83) 
The mean value (Ha) 2 is then given by 

00 

(H,), = (B„ H^B,) = ^^(BjEfc = \z\^ + 1 (84) 

fc=0 

where we have used Eq. ( 1751) . In the same manner we find (H^)^ = |-2|^ + 4|zp + l, so that 
AHa = \z\\pl. Hence, for very large \z\ one gets AH^ << (Ha)^ and the relative value of 
the energy of the state B^ is well defined, as usual for the Glauber states. It is also simple 
to verify that the product of the root-mean-square deviations AP and Ay is minimized. 
In conclusion, the 2;-parameterized functions (182|) are the coherent states belonging to a 
wide class of position-dependent mass oscillators of the first kind. 
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5 Concluding remarks 



We have studied the problem of solving the Schrodinger equation for an arbitrary position- 
dependent mass system. Our approach is useful to face two general physical situations. In 
the first one we look for the interaction which must be applied on a mass m{x) to supply 
it with a particular spectrum of energies. The second physical situation corresponds to 
the case in which one is interested on a given position-dependent mass m{x), subject to 
a particular interaction V{x) rather than in the recovering of a specific spectrum. For 
arbitrary orderings of m{x) and P in the Hamiltonian, diverse general expressions for 
m{x) were derived as a consequence of mapping the original Schrodinger equation to a 
conventional constant mass one. It was also found that the transformation is rather simple 
for the very special ordering m~^^^P'm~^/'^P'm~^^^ in the kinetic part of the Hamiltonian. 
In contradistinction with [21,22,28], we showed that a position-dependent mass Hamilto- 
nian can be factorized as the product of two mutually adjoint operators with no a priori 
assumptions on the ordering of m and P. 

In particular, two kinds of position-dependent mass oscillators were analyzed. The 
first one is defined to be isospectral with the quantum oscillator of mass mo, no matter 
the explicit form of m{x) and V{x). The oscillators of the second kind exhibit spectra 
different from the equidistant energies hwo{n + 1/2) and correspond to a particular mass 
m(x) subjected to the harmonic oscillator potential. Results include the singular oscillator 
as well as confining odd-root-law, In^ and sinh^-like interactions. The special case of a 
particle of mass tjiq in a confining even-power-law potential was found to be isospectral 
to a system of mass m{x) oc x"^"/*^^""*"^), subject to the action of the harmonic oscillator 
potential. The factorization operators were then selected to work as ladders in the space 
of the position-dependent mass oscillators. Finally, the coherent states corresponding to 
oscillators of the first kind were explicitly constructed as eigenvectors of the annihilation 
operator. These new CS have the same analytical form as the Glauber states and minimize 
the position-momentum uncertainty principle. 

Of special interest, the singular oscillator V^gq defined in Eq. (1481) and referred in 
Section [3.2.3l as the squeezed oscillator, exhibits CS connected with the su{l, 1) Lie algebra 
if the mass is a constant [30]. Particular cases of the mass-function have been shown to 
preserve the su{l, 1) spectrum-structure of V"sq [28]. The same is true for any of the masses 
derived in this paper. Thereby, it is sound to construct position-dependent mass su{l,l)- 
hke CS. Following [30], such result could be applied to get a better understanding of the 
physics of N interacting particles (work in this direction will be published elsewhere [48]). 
Other physically interesting systems can be analyzed in the corresponding manner once 
the dynamical algebra is given. Special attention must be drawn to the Susy non-linear 
algebras engaged with infinite point spectra. If the energy levels can be obtained by a 
function of its index En = E{n) then one can distinguish between natural and linear 
algebras of the Susy-system. To each one of these algebras there exists a companion set 
of CS [32]. Then, besides the systems above discussed, it would be also interesting to 
analyze the position-dependent mass CS belonging to higher-order Susy potentials like 
the Poschl- Teller ones [49] and [50] (see also [36]). Results on the matter are in progress. 
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